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Fermi, Pasta and Ulam observed, that the excitation of a low frequency normal mode in a nonlinear 
acoustic chain leads to localization in normal mode space on large time scales. Fast equipartition 
(and thus complete derealization) in the Fermi-Pasta-Ulam chain is restored if relevant intensive 
control parameters exceed certain threshold values. We compare recent results on periodic orbits 
(in the localization regime) and resonant normal forms (in a weak derealization regime), and relate 
them to various resonance overlap criteria. We show that the approaches quantitatively agree in 
their estimate of the localization-delocalization threshold. A key ingredient for this transition are 
resonances of overtones. 
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I. INTRODUCTION 



Fermi, Pasta and Ulam (FPU) considered the dynamics of a classical nonlinear acoustic chain of particles with 
coordinates x n , canonically conjugated momenta p n and the Hamiltonian 
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where a and (3 are the nonlinearity parameters. Fixed boundary conditions x Q — sjv+i = Pa = Pn+1 = are assumed. 
N is the number of degrees of freedom of the system. FPU considered the case with a > and (3 = 0, and the one 
with a = and (3 > 0. These two models fll} are known as the a-model and /3-model, respectively. 

FPU performed a numerical experiment [l[ with the aim to observe an approach to an expected energy equipartition 
among the Fourier modes of system ([T]), when the longest wavelength mode is initially excited. Instead, they discovered 
that the energy, initially given to the longest wavelength mode, is partially transferred just to a few modes of shorter 
wavelength that are excited one after the other in order of decreasing wavelength. Such an energy cascade was 
observed to be effective, in practice, up to a certain critical mode, giving thus rise to a state of the system quite 
different from the equipartitioned one and apparently persisting up to their largest integration time, in the range of 
parameters (energy, number of degrees of freedom, nonlinearity parameters) they chose. Such an observed lack of 
complete energy equipartition among the modes of the system up to unexpectedly long times posed a problem which, 
after them, was named the FPU problem or paradox. We stress here, that the observation of an energy distribution 
in normal mode space staying localized for just the times on which nonequipartition is registered, is perhaps the main 
issue of the paradox. Indeed long-studied questions of estimating recurrence times are already based on the fact that 
the energy stays localized in normal mode space for long times Q. Galgani and Scotti Q observed numerically in 
1972, that the localization profile is exponential. 

The basic mechanism of the energy cascade in the FPU problem was identified by Ford (1961) [4(. The long 
wavelength modes are characterized by a dispersion relation of the type uj q ~ nq/N, where q = 1, . . . , N is the mode 
index, so that 



LOl ~ L0 2 /2 ~ W 3 /3 . 



(2) 



Due to such an almost complete acoustic resonance condition a resonant excitation of consecutive modes of decreasing 
wavelength, or increasing frequency, is expected. Notice that the approximate equalities ([2]) hold with an error 
increasing with q, so that one can expect the resonant pumping mechanism to become uneffective at a certain critical 

In 1966 Izrailev and Chirikov (IC) [5] realized that increasing the energy of the initially excited mode will lead to a 
crossover from the localization (nonequipartition) regime to a delocalization (equipartition) one, where a fast transfer 
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of energy from long to short wavelength modes will lead to energy equipartition, as originally expected by FPU. IC 
used a nonlinear resonance overlap criterion to estimate the crossover. Notably they considered resonances different 
from (|2l). 

Both the work of IC and the idea of Ford have been re-considered by Shepelyansky [6( in 1997, who applied the IC 
idea by using the Ford resonances ([2]). He obtained a resonant normal form for the system {T]), from which he arrived 
at an estimate of the critical mode index q c at which the cascade effectively slows down. Numerical runs yielded 
that the numerical distribution of the harmonic energies E q of the modes (i.e. E q vs. q) displays an exponentially 
decreasing tail of the form E q ~ exp(— q/q c ). Since Shepelyansky used resonances which are different from the ones 
originally used by Izrailev and Chirikov, the crossover estimate from [fjj does not coincide with the IC result. 

Two more recent lines of research are the reason for this report. First, a normal mode at the linear limit a — [3 — 0, 
when excited to a given energy, represents a periodic orbit in the phase space of the system ([T]). That periodic orbit 
persists for nonzero a, B at the same energy and is coined g-breathcr. Existence proofs, perturbation theory and high 
precision numerical calculations @, H, H, Hcj, EL 03 allow to draw a consistent picture of these periodic orbits. For weak 
nonlinearity q-breathers are exponentially localized in normal mode space. Analytical estimates allow to conclude 
about the existence of a derealization threshold in the control parameters, above which the (/-breather delocalizes. 
Small deviations from these orbits stay close to them for sufficiently long times. Thus some features of the trajectory 
first computed by FPU are quantitatively captured by the properties of a g-breather and by the linearized phase space 
flow around it. The scaling properties of g-breathers appear to be close to the results of Shepelyansky. 

Second, a careful reconsideration of the resonant normal form [l3l [Til [TEl. ITrl ITtI [l8j allowed to get analytic estimates 
of the modal energy spectrum of the system and of the time scales involved in the problem. The direct relevance of 
the work of Zabusky and Kruskal on the KdV equation to the FPU problem has been almost completely clarified. 
Again some of the scaling properties of the normal form are similar to those of Shepelyansky, and thus related to 
the scaling of q-breathers. That is remarkable, since q-breathers and normal forms seem to operate in conceptually 
different regimes, namely the localization and weak derealization regime, respectively. 

Here we intend to compare all these results with each other, in order to establish the similarities and differences. 
The emerging picture is, that the relaxation properties of a normal mode, launched into a nonlinear acoustic lattice, 
are governed by intensive quantities (i.e. energy density, frequency, etc). Below certain parameter thresholds the 
relaxation is slow, while above them it is fast, and that circumstance may survive in the thermodynamic limit. 



II. TRANSFORMATION TO NORMAL MODES 



The usual canonical change of variables (x,p) i— ► (Q,P) that transforms Hamiltonian (JTJ) into the sum of a Hamil- 
tonian of nonintcracting mode oscillators plus a perturbation responsible for mode coupling is given by 



where 
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for q = 1, . . . , N. In terms of the canonical normal mode variables (Q, P) the Hamiltonian ([TJ) reads 

N P 2 + LU 2 2 N 

H(Q,P) = J2 Q 2 + 3V27V + 2 ^ S 3 (j>Q,l)^^lQjQ q Ql + 

9=1 i, 

3 N 

+ 8(N + 1) ^ S i(^^^ m ) U i iJ <l iJ l U> rnQjQqQlQm , (5) 
j,q,l,m— 1 



^ =2sin (^y (e) 

is the dispersion relation (frequency vs. mode index) of the normal mode of the noninteracting system, whereas the 
two coefficients 5*3 and S4 appearing in ([5]) rule the actual mode coupling in the system, and are given by 

Sz(j, q, I) = 8j+q,l + 8j+l,q + Sq+lJ - Sj+q+l,2N+2 , (7) 
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^(j, 9, I, ™) = <5j+ij+f,m + 5 q+ i +m j + 5l+ m +j t q + 5 m +j+q,l + 

Sj+q.l+m + <5j+2,q+m + <5j+m,<}+2 — &j+q+l+m, 2N +2 + 

~<5j+9+2,m+2JV+2 — ^q+l+m,j+2N+2 ~ 3l+m+j,q+2N+2 — <5m+j+g,i+2JV+2 > (8) 

<5n,m denoting the usual Kronecker symbol. The coefficients S3 and 5*4 are invariant under permutation of any pair of 
indices. Notice that the dispersion relation ©, when q/(N + 1) is small, admits the expansion 
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0((q/Nf) , (9) 



so that, in a very first approximation, low frequency modes (small q) display a quasi-linear dependence of the frequency 
on the mode index and their frequencies approximately satisfy the relation ([2]). The equations of motion associated 
to the Hamiltonian ([5]) read 
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Qq+^lQq = --y==^2 S 3 (q,j,l)u} j (JiQ j Qi + 
/3uj N 

~2N+~2 ^ S -±(<l'j> l >' m ) UJ 3 UJ l UJ ?nQjQlQ,n ■ (10) 
j,l,m— 1 

They are used for constructing g-breathers. For what concerns the normal form approach, a further canonical transfor- 
mation to complex variables turns out to be useful. More precisely, one defines the change of variables (Q, P) i— > (z, z*) 
by 

Zq = + (11) 

W2u) a 



and its complex conjugate z* for any q = 1, . . . , N. Notice that 

E q ^u q \z q \ 2 = {P 2 +uj 2 q Q 2 q )/2 (12) 

is the harmonic energy of mode q, so that \z q \ 2 is the action variable associated to the same mode. Regarding z q as 
a canonical coordinate, it turns out that its conjugate momentum is given by iz*, since 

{•^rr %ti } — i^r.s , 

where {...,...} denotes the usual canonical Poisson bracket. Thus the Hamiltonian ([5]) becomes 

JV N 
H = '^2 UJ q \z q \ 2 + — y= = 5 ' 3 ^' 9 ' l)V^jUq~Ul{Zj + Z j)( Z q + Z q)( Z l + Z l) + 

9=1 i,9,i=l 
£ N 

+ 32(N + 1) ^ Si<y ^ q ' l ' m )V^W^( Z J + Z J )( Z 1 + Z q)( Z l + + z *n) > (13) 

j.q,l,m—l 

with equations of motion given by iz q = dH/dz*. 



III. THE FPU PROBLEM 



Fermi, Pasta and Ulam launched a trajectory by initially exciting just one normal mode with seed mode number 
go, and monitored the time evolution of the harmonic mode energies (|12p. A typical outcome of their experiment is 
shown in Fig[TJ Instead of spreading among all normal modes, the energy stays localized on a few long wavelength 
modes. Recurrence events are observed on times much shorter than the expected Poincare recurrence times. At the 
same time the localization of energy in normal mode space happens for times which are orders of magnitude larger 
than the recurrence times. 

In Figf5]we show distributions of normal mode energy densities e q — E q /N for larger initial energies, and at different 
times t = 10 4 , 10 5 , 10 6 . To eliminate fluctuations, they are averaged over a preceding time window of 10 4 . We clearly 
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FIG. 1: Normal mode energy evolution for q = 1,2,3,4,5 for (i) an FPU trajectory (oscillating curves) and (ii) an exact 
q- breather periodic orbit (horizontal lines). Parameters are a = 0.25, E = 0.077, N = 32 and q = 1. From @] 
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FIG. 2: (a) Distribution of normal mode energy densities for an FPU trajectory and parameters qo = 1, N = 31, a = 0.33, 
E — 0.32. Circles - t = 10 4 , squares - t = 10 5 , rhombs - t = 10 6 . Dashed line - distribution for the g-breather periodic orbit 
from (b) for comparison, (b) Distribution of normal mode energy densities for a g-breather periodic orbit and the parameters 
as in (a). From hj 
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observe exponential localization of the distribution, and note that such a distribution is reached on a quite fast time 
scale t \ <C 10 4 . We also observe a slow resonant excitation of tail modes, which ultimately drags the system to 
equipartition, but on a much larger second time scale Ti ^> T\ [2(| ■ The origin of these tail resonances is explained in 

Further increasing e.g. the seed mode energy, or the nonlinearity, may lead to a threshold where fast derealization 
and equipartition take place, i.e. where ri ~ ti [2(3] • Below these thresholds, the two time scales are strongly 
separated, and the evolution on times ri < t < n corresponds to the regime observed by FPU. Such a double 
time-scale scenario is akin to that proposed first in plj . From this perspective, the FPU problem in a nutshell may 
be phrased as: Why are there (at least) two different time scales, how to estimate the localization length and what 
is its relation to the time scales, is the effect robust in the limit of infinitely large system sizes, as well as for higher 
spatial dimensions? 



IV. CONTINUATION OF NORMAL MODES: g-BREATHERS 

Let us assume, that the FPU trajectory is close to a periodic orbit (PO) for times i«T 2 . Then the FPU trajectory 
evolution will be almost regular, as observed in the simulations. It can be expected, that the PO will correspond to 
a slightly deformed normal mode solution of the linear case [13] • 



A. General properties 

A rigorous construction scheme of such POs (including an existence proof) was given in [7, 8], which starts from the 
simple linear limit a = (3 = 0. Essentially one chooses a seed mode number go and an energy E. This choice uniquely 
defines a periodic orbit in the phase space of (U) for a — (3 = 0. Moreover, for any finite size N the nonresonance 
condition nu> qo ^ ow Qn holds for any integer n [8j . Therefore straight application of a persistence theorem of periodic 
orbits by Lyapunov [231 ] guarantees the continuation of that periodic orbit for nonzero a and j3. Notice that for 
uniqueness one has to fix one parameter, e.g. the total energy E. 

The continued periodic orbits are exponentially localized in normal mode space and follow the corresponding FPU 
trajectories impressively closely (cf. Fig [1] and Fig They are coined g-breathers @> H B E3, EH E3] • The localization 
length tends to increase with increasing nonlinearity, or increasing energy, or decreasing seed mode number. Numerical 
high- precision schemes for computing these orbits are given in @, Q . 

Notably g-breathers may stay stable, or also turn unstable, keeping their localization 0, Q. Such instabilities, 
which are related to the emergence of weak chaos in the core of the corresponding FPU trajectory [24| . do not affect 
the time scale T2 to a significant extend. That is reasonable, since an instability marks a resonance, but the fact that 
the g-breather is well localized tells then, that the resonance can not be between core and tail (which would lead 
to derealization in the first place). Thus the only possibility is a resonance between two modes in the core, which 
changes the core dynamics, but leaves the localization properties almost unchanged. 

The concept of g-breathers has been successfully applied to higher dimensional systems as well, without much 
change in their properties @. We think that most of the results to be discussed below, though strictly speaking being 
correct only for one spatial dimension, will hold in higher dimensions as well. 



B. Scaling to large sizes 

Application of Poincare-Lindstedt perturbation theory to g-breathers 0, H, H[ suggests that only intensive quantities 
matter, namely the energy density e and the wave number k defined by 

e = E/N , k = nq/(N + 1) (14) 

matter. It follows that g-breathers may persist in infinitely large systems. As was shown in [Io| . it is enough to first 
obtain a g-breather for a chain with some finite N. Let us then consider a new chain with size N + 1 = r(N + 1), 
with r = 2,3,4, .... It follows, that 

Qdt) = {f QM (15) 

is a solution of the scaled chain. Fixing r and repeating that procedure, or simply considering the limit r — > oo, we 
will obtain solutions for macroscopic systems. These scaling procedures are easily generalized to two and three spatial 
dimensions, as well to other boundary conditions [lfj| . 
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FIG. 3: S m {z) (dashed line). Symbols and guiding line: data from (/-breather calculations for systems with /3 = 1, e = 
6.08 ■ 1(T 4 , 9.6 ■ 10" 4 , 1.57 ■ 10" 3 and N = 149, 359. From [l(J 



C. Localization properties: j3 FPU 

The perturbation theory results @, [§[ for the FPU /3-model yield the following form for the modal energy spectrum 
of a q-breather with seed wave number fco [To| : 



In e k = - 1^ In VA + In e ko , VA = 



8 kl 



(16) 



Here e ko is the seed mode (harmonic) specific energy. Localization holds for A < 1 , which is a function of intensive 
quantities only. It takes an exponential law 

e k ~ e- fe /« (17) 



where £ is the localization length (in fc-space). We consider g-breathers at fixed average energy density e. Using (JTJ 
it follows ek = (1 — A)e. Together with the definition of A in (TlTJl) we find, that the inverse of the localization length 
£ in fc-space is given by: 

T^T^nVA, VA=^(V / TTIF^-1), z=^, - 2 = ^e, (18) 

where z is a scaled wavenumber and v the effective nonlinearity parameter. It follows, that £ _1 = —S m /v where the 
master function S m = z^ 1 lnVA depends only on the scaled wavenumber z. S m vanishes for z — > and has its largest 
absolute value max(\S m \) ~ 0.7432 at z m i n m 2.577. For a fixed effective nonlinearity parameter v the g-breather with 
k m in = vz m i n shows the strongest localization, while larger and smaller fco tend to weaken the localization. Especially 
for ko — > the (/-breather delocalizes completely. With increasing v : the localization length of the g-breather for 
ko = kmin increases. For k Q ^> v it follows S m « — 2/zln(z) and for k <C v we find S m » — z (see FigEJ). All 
predictions are verified numerically. It follows, that in the regime of strong localization e w e ka the exponential 
decay in units of fco depends on the single parameter 3/3e/c /(8fco), implying corresponding scaling laws. In particular, 
delocalization is expected to be reached when 



However, for long wavelength it follows 



3/3 e 

P 
o 



fir 1 - (19) 



£ = v 2 /ko , A; < v . (20) 



Since fc-space is finite, localization is only meaningful if (at most) £ < n. Consequently g-breathers delocalize at a 
threshold given by 



Clearly the scaling law is now different from the one for strongly localized ^-breathers (|19p . 
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The reason for the weaker localization of g-breathers when ko ^> v is the increasing distance 3fco between modes 
excited in consecutive orders of perturbation theory. The delocalization for fc however is due to an approaching 
of resonances (2n + l)wfe — > W( 2Tl +i)fe for some integer n, which are the Ford resonances j2]). 

Most importantly, for each value of (3 and e we can identify a critical seed wave number 

kd = ^ ■ (22) 

For fco > feci g-breathers are localized, while for fco < fc c i they delocalize. This delocalization threshold, corresponding 
to condition (f2"Tjl . can be expected to translate into a delocalization threshold of a corresponding FPU trajectory. If 
instead we use the scaling law (|19p for strongly localized g-breathers, we arrive at a critical seed wave number 



¥ • (23) 



D. Localization properties: a FPU 

In a similar study [12] the a-FPU case yields analogous results. The modal energy spectrum of a g-breather with 
seed wave nuber ko is given by 



1 (£- i )" ,7++21 °(£)""'"' -<- aj w- 



Note that the exponential decay here is dressed by a weaker but nonzero power law, which may be important when 
it comes to quantitative comparisons. Using (|24[) it follows efe = ^"^7 e and thus 

-4( 1- 7) 3 k r 1/4 

7 = z — — , z = — ,(/> = y/a e 1 . (25) 

1+7 

Here we again introduced a scaled wave number z and the effective nonlinearity parameter (f>. The inverse localization 
length £ of a g-breather is then defined by a master slope 5 m via 

r X = -%, S m = iln 7 . (26) 

The dependence of S m (z) is qualitatively similar to the /3-FPU case. It has a minimum at z m i n « 2.39 and a value 
S m (z mm ) w —1.5. At variance with the /3-FPU case the dependence for small z <C 1 is S m ~ — z 1//3 . 

All predictions are again verified numerically. It follows, that in the regime of strong localization e w £fc the 
exponential decay in units of fco depends on the single parameter a 2 ek /kQ, implying corresponding scaling laws. In 
particular, delocalization is expected to be reached when 

^ = 1 . (27) 

However, for long wavelengths it follows 

£ = 4/3 /(2fco) 1/3 , k « <f> . (28) 

Since fc-space is finite, localization is only meaningful if (at most) £ < tt. Consequently g-breathers delocalize at a 
threshold given by 

2fc 

Again the scaling law is now different from the one for strongly localized g-breathers (|27p. 

The reason for the weaker localization of g-breathers when ko ^> <fi is the increasing distance 2ko between modes 
excited in consecutive orders of perturbation theory. The delocalization for ko — > however is due to an approaching 
of resonances nuk — > co n k for some integer n > 2, which are the Ford resonances |2|). 
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Most importantly, again for each value of a and e we can identify a critical seed wave number 



a 2 e 



fc* = ^3 • (30) 



For fco > fed g-breathers are localized, while for fco < fcci they delocalize. This derealization threshold (|2 ljl is expected 
to translate into a derealization threshold of a corresponding FPU trajectory. If instead we use the estimate (|27[) . 
we arrive at a critical seed wave number 

k c2 = aV2 e l/4 . ( 31) 



E. Some consequencies 



Let us discuss the connection between the time scales n and t 2 , and the localization-delocalization threshold for 
g-breathers. According to our understanding, for fc > fc cl it follows t 2 3> T\. For fc « k c i we reach r 2 ~ T\, and 
when further decreasing ko < k cl , both scales stay close, but diverge in the limit of long wavelength. Note that 
the critical number k c \ has to be replaced with k C 2 if one takes the scaling results for strongly localized g-breathers. 
These conclusions are only reasonable if we start a trajectory close to the g-breather. Thus, if the g-breather is almost 
(or completely) dclocalized, an initial condition similar to the one considered by FPU (which is strongly localized in 
g-space) is at large distance from the corresponding g-brcather, and the relaxation properties of the FPu trajectory 
may become very different. 

For a large (possibly infinite) system, at some given energy density and nonlinearity parameter, a g-breather 
necessarily delocalizes for some small but finite seed wave number fc , signalling the breakdown of the corresponding- 
dressed phonon picture familiar from solid state physics. That may seem strange, since it is textbook wisdom that 
acoustic phonons have lifetimes which increase with increasing wavelength. Yet the two results perfectly match. 
Indeed, the lifetime measures the relaxation time of a bare phonon (normal mode) injected into the system, in units 
of its period. These lifetimes indeed grow with increasing wavelength, since the coupling to shorter wavelength modes 
is inversely proportional to the third (a) or fourth (/3) power of the wavelength. At the same time the couplings 
become more and more resonant, so that the new state which is reached after that waiting time becomes more and 
more different from the bare phonon state (normal mode). 

Finally, we emphasize that quantum phonon physics usually considers the decay of one phonon into two phonons 
of longer wavelength. However in the classical limit, considered here, the occupation number of the initial phonon 
state is large compared to one, and the leading order decay is given by annihilating two (a) or three (/3) phonons and 
creating one with corresponding shorter wavelength, in full agreement with the numerical simulations of FPU, and 
with all the other studies of the corresponding classical problem. 



V. THE NORMAL FORM APPROACH TO THE ENERGY CASCADE 



The picture emerging from the above mentioned approach to the FPU problem via g-breathers leads to the conclusion 
that the perturbative continuation of a normal mode of wavenumber kg of the linearized system ([T]) converges to an 
exact periodic solution of the full nonlinear system when fco is larger than some boundary or critical value fc c , the 
latter quantity being a function of the specific energy e of the system. Under the same condition, fco > fc c , the modal 
energy spectrum of such a periodic orbit is exponentially localized in fc-space. Thus, in the corresponding initial value 
problem, when a mode of wavenumber fco > fc c is initially excited, the FPU paradox can be explained in terms of 
closeness of the initial datum to an exact periodic orbit of the system. On the other hand, the continued orbit of a 
mode with wavenumber fco smaller or much smaller than fc c , when it exists, is expected to be delocalized in fc-space, 
and the dynamics of the corresponding initial value problem, with initial excitation of a mode with wavenumber 
fco <C fc c , becomes different. In the latter case the energy initially given to the mode with fco <C k c is expected to be 
transferred to higher wavenumber modes within some characteristic time r c . 

The qualitative outcome is sketched in Figf4j There exists a wavenumber interval in fc-space, tt/N < fc < k c (e), 
that one can refer to as the acoustic resonance range. The energy initially injected inside such a wavenumber range 
cascades up to the border fc w fc c , within a characteristic time t c , at least if the specific energy e is small enough. The 
subsequent transfer of energy inside the range fc c (e) < fc < 7r, that one can refer to as the dispersive range, is slowed 
down by the detuning from the almost complete acoustic resonance of the corresponding unperturbed frequencies. 
The path to energy equipartition is expected to be completed on some time scale(s) much longer than r c . Such a 
picture emerged for the first time in [y, [25l |26| . where estimates of fc c (e) and of T c (e), were provided, mainly for 
the a-model. Moreover, it was numerically shown there that the modal energy spectrum of the system displays a 
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FIG. 4: A synoptic diagram describing the energy cascade in FPU (the thick solid curve is the graph of the dispersion relation 
©). The vertical thin-dotted line marks the position of the critical wavenumber k c . The long-dashed straight line is the graph 
of the dispersion relation linearized around k = 0, which yields the Ford resonances ([2]). The fast energy cascade takes place 
on the time scale r c (e), whereas a second slower stage of the cascade is expected to continue into the more distant part of the 
wavenumber space, as observed numerically. 
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Wien-like exponentially decreasing tail inside the dispersive range, namely E(k) ~ e ~k/k c f or ^ ^ ^ ag originally 
observed in [3J and stressed again in [2l| . 

The fundamental scaling laws of the problem (i.e. the dependence of k c and r c on e) can be obtained through the 
approach based on the quasi-resonant normal form construction introduced in [|| and developed in [3, EE EE E3, El] • 
The idea starts from the expansion ^ of the dispersion relation ©, namely cu q — 2sin(fc 9 /2) = k q — fcg/24 + . . . 
(recall here that k q = irq/(N + 1), where q — 1,2...), which holds if k q is small. The quadratic part hi of the 
Hamiltonian (fl3|) can thus be rewritten as h 2 = Ylq=i LJ q\ z q\ 2 = J2 q =i kq\ z q\ 2 + S^Li ^ 9 1 -^q 1 2 , with Q q = uj q — k q . For 

initial excitations of acoustic modes, the sum Eh = Y^q=i kq\ z q\ 2 can be regarded as the unperturbed Hamiltonian 
of the problem. Such a splitting holds if low modes only are involved in the dynamics and consistency is recovered 
at the end, showing that the fraction of modes k c sharing the energy turns out to be small when the specific energy 
e is small. One now introduces a time-dependent canonical change of variables z h- > £, the new variables £ q (t) being 
defined by z q — C, q exp(—ik q t) (g = 1, . . . , 2V), in terms of which the Hamiltonian of the problem is mapped to a new 
Hamiltonian K = K((, (*) + R(C, (*,t), where 

N N 



9=1 j,9,i=l 

8(JV + 1) 



+ 8(N+ 1) ^ WW^^" (CiC0*Cm + c.c.) 4- 
j.q,l,m=l 



3 p N 

^ Sj+q,l+m^UJjUJqUJlUJ m (CjCqClCm + C.C.) (32) 



32(JV+1) . , 

j , q , / , m — 1 
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and the remainder R collects all the terms that explicitly depend on time through some oscillating exponential. The 
time-average R of R, the Cfc's being held constant, is zero, and one can study, in a first approximation, the dynamics 
associated to the Hamiltonian (f32f which is precisely the quasi-resonant normal form Hamiltonian of the FPU problem. 
The rigorous justification of such a n ap proximation, i.e. the neglecting of the remainder i?, is a difficult task, and 
precise statements are found in [l?], Il8j |. Before passing to a direct inspection of the equations of motion associated 
to the Hamiltonian (f32|) for the a-model and for the /3-model, we observe that, as a consequence of averaging, the 
unperturbed Hamiltonian Eh = „_ 1 fcg|Cg| 2 commutes with the normal form Hamiltonian (|32p . The second integral 
Eh approximately coincides with the sum of the harmonic energies of the modes, at least for excitations involving 
acoustic (low frequency) modes only. 

Let us now consider the a and (3 models separately. In both cases, one supposes that only a fraction of the order 
k c <C 1 of acoustic modes is involved in the dynamics, so that ( q ~ for q > Nk c (k c to be estimated below). One 
then introduces the new noncanonical variables 



in terms of which the second integral Eh, i.e. the sum of the acoustic harmonic energies, becomes 



N 

E h = ^2 \u k 
fc=i 



2 _ 



Ns 



(33) 



(34) 



In the latter formula we have defined the harmonic specific energy e as the acoustic harmonic energy per degree of 
freedom. Of course, as far as the perturbative approach holds, one has e ~ e. Notice that the quantity |u g | 2 is the 
approximate harmonic energy of the mode q. 



a— model 



The equations of motion of the a-model, in terms of the u variables (|33[) take the form 



k 3 
fi 

24 



k q a 



4ViV+T 



9-1 

U=l 



N-q 

1=1 



The fundamental step consists now in recognizing that, through the scaling transformation (u q ,k q ,t) 
defined by 

u q =e l ' 2 v q , k q = c^V/ 4 ^ , t = a~ 3 / 2 e^ 4 T , 
the equations of motion (|3"5|) and the second integral (|3~4"1) transform to 



.dv q 



24 



q-1 N-q 
^Vq-lVl +2^2 V*lV q+ l 
.1=1 1=1 



1 N 

£2> 



(35) 

[Vq, X q , T"), 

(36) 

(37) 
(38) 



8=1 



The latter system no longer depends parametrically on a and/or e, the only parameter entering the equations of 
motion being N. We now define an effective Reynolds number lZ(x q ) as the ratio of the nonlinear to the linear 
term appearing on the right hand side of (|37[) . both terms being squared and averaged over the phases of the active 
modes (see (14j). Thus TZ(x q ) measures, for any mode q, the relevance of mode coupling (nonlinearity) with respect 
to dispersion. It follows that, under the hypothesis that the phases of the active modes are randomly distributed, lZ(x) 
decreases below one at x c ~ 1, independent of N. From the equations of motion (|3T[) it follows that the time scale 
T c associated with the dynamics of the mode with rescaled wavenumber x c is given by the reciprocal of its frequency, 
namely T c w (x c )~ 3 . Thus also T c ~ 1 does not depend on N and, from the second and third of the definitions 
the estimates 



\fa e 



1/4 



(39) 
(40) 



follow. Notice that the scaling of k c coincides with that of k C 2, equation (f3"Tj) . 
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B. /3-model 



For what concerns the /3-model we will write down the main formulas, since the reasonings are much the same as 
those made for the a-model. The equations of motion of the system written in terms of the u variables (|33[) are 



kg 

24 



k q (3 



8(JV + 1) 



^ UlUmUq-l- m + 
2<l+m<q-l 



+3 ^2 u *lKi U l+™+q +3 U* l+m _ q UlU m 

2<l+m<N-q q+l<l+m<N+q 



(41) 



while the conservation law (|34[) is exactly the same. The scaling transformation (itk ,£,k,t) ^ (vk ,Xk,T) is now given 
by 



,1/2 



v q , £ 9 = (/fe) 1/2 x q , t=(/3e)- 3 / 2 T , 



(42) 



through which the equations (|41[) transform to 



.dv q 



Z-3 

24 



Vq+ 8(N + 1) 



2<l+m<q-l 



vfv^Vl+m+q + 3 Y 



l+m — q 



VlV r , 



2<l+m<N-q 



togheter with the conservation law 5Z g=1 \v q \ 2 /N = 1. The ratio nonlinearity/dispersion lZ(x q ) is then defined in a 
way completely analogous to what done above for the a-model and, under the hypothesis that the phases of the active 
modes are randomly distributed one gets that IZ(x) < 1 as x > 1. Thus, also for the /3-model x c ~ 1, T c ~ l/(x c ) 3 ~ 1, 
both independent of N, and from the scaling transformations (I42|) one gets the estimates 



q+l<l+m<N+q 



(43) 



t c ~(/3 £ )- 3 / 2 . 

Here again we stress that fc c displays the same scaling of fc C 2 given in ([23 



(44) 
(45) 



Comments 



First of all, note that the obtained scaling properties of k c correspond to the scaling properties of k c i for ^-breathers, 
for the a and (3 models, respectively. This is quite natural: the region where g-breather miss their strong localization 
property is the one where the acoustic resonance effectively works, leading to a partial relaxation of the system. Recall 
however that k C 2 is obtained using the scaling relationships for strongly localized g-breathers, while a more careful 
analysis yields estimate k c i which scale differently and are lower than k C 2 . That leaves room for further investigations 
and comparisons. 

As a further comment, making reference to Figf4l we stress that in order to observe the scenario depicted above 
(fast cascade up to k c and slow one from k c on) it is necessary that ir/N <C k c (e) <C 7T, so that the acoustic resonance 
range has a small but finite size, and the dispersive range occupies the most part of fc-space. The upper bound on k c 
gives a threshold value of the specific energy above which a fast trend to equipartition occurs. The lower bound on 
k c yields instead an estimate of the minimum number N of degrees of freedom necessary to open an acoustic range in 
the system at a fixed total energy E. For the a-model one gets N 3> q; -2 / 3 ^" 1 / 3 , whereas for the /3-model it must 
be N » l/(/3E) (neglecting numerical factors). 

We want to stress here that the hypothesis of random phases seems to play a crucial role in the above derivation, 
in particular for what concerns the dependence of k c and r c on the number N of degrees of freedom of the system. 
We anticipate that when one replaces such an hypothesis with that of coherent phases, for example, the estimates 
([3"5)) - (|4"5"l) hold the same but for the replacement of the specific energy e with the harmonic energy Eh — Ne, which 
approximately coincides with the total energy of the system. Whether the phases of the active modes are randomly 
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distributed or not seems to depend strongly on what one chooses initially: if one initially excites a small fraction of 
acoustic modes, randomness or coherence in the phases of the modes initially excited is recalled by the system, and 
such a strong memory effect is another signal of quasi- integrable behaviour [271 ]. 

It must be stressed that the scaling (|39p . first obtained by Shepelyansky, both numerically and analytically, has 
been confirmed numerically in two successive works by Berchialla et al. [25[ and Biello et al. [261 ] . Notice that the 
characteristic time t c (|40j) was first written down in [6| on the basis of dimensional arguments and interpreted as 
the largest Lyapunov exponent of the system, whereas it was first measured in [26| and identified as the time scale 
at which energy is transferred to short wavelengths. Moreover, the analytical explanation of the exponential tail in 
the modal energy spectrum has been given in [l3|, EH EE El, E3, EB] , where it has been shown that the FPU quasi- 
resonant normal form dynamics in real space is described by two KdV equations for the a-model and by two modified 
KdV equations for the /3-model. Peculiar integrability properties of the KdV equations yield then an exponentially 
decreasing modal energy spectrum; see in particular [13j where, based on soliton theory, a formula describing all the 
spectrum of the a-model (not only the high frequency tail) is given. 

Another consequence of the analysis carried out in terms of normal forms is the following. At least in the case of 
the a and /3-models, the normal form being integrable, it turns out that the first stage of the dynamics, namely up 
to times of the order of r c , is almost regular. Thus r c has to be interpreted just as a redistribution time, and not as 
a stochastization time scale. As a good example one can think of the Toda model, whose dynamics, for low specific 
energies, displays all the features of the a-model. However, such a similarity breaks down on some time scale, after 
which the dynamics of the Toda chain goes on to be regular (quasi periodic with zero Lyapunov exponents) while 
that of the a-model develops chaotic features (local stochasticity and a nonzero Lyapunov spectrum); see [28|. The 
onset of stochasticity in the FPU dynamics takes place on time scales surely longer than r c (e). 



VI. RESONANCE OVERLAP CRITERIA 



The idea of the resonance overlap is to obtain parameter estimates where a certain resonance in the nonlinear 
system may develop, leading to a strong change in the dynamics of the chain, presumably to a fast equipartition. The 
main ideas are published in [f| and (29| , and we are going to present here a very simple and qualitative way to derive 
the final results for the /3-FPU case. 

The selfinteraction of an initially excited mode with seed wave number kg is described by 

Ofco+<~f<Qt - (46) 

This is an anharmonic oscillator problem, where the frequency depends on the energy (30j . Consequently the frequency 
shift of this anharmonic oscillator problem, when excited to a given energy E — Ne, is given by 

8uj~f3uj ko e. (47) 

When this shifted frequency matches some other frequency, a resonance may develop, leading to a resonant decay of 
the originally excited mode. The entire final result will now depend on the choice of the frequency to be compared with, 
which depends on the model. Izrailev and Chirikov chose the seemingly natural frequency given by the separation of 
the normal mode frequencies due to the finite size of the system. Such a frequency separation is Ajc ~ 1/iV for long 
wavelength modes. Comparing it with the nonlinear frequency shift 5oj one arrives at the threshold condition 

fak e ~ i ■ (48) 

Clearly this condition is not entirely depending on intensive quantities. Furthermore no correspondence to the above 
derived derealization thresholds for g-breathers and scaling relations for normal forms can be found. Finally note, that 
for the longest wavelength mode q = 1 it follows that the critical energy density according to IC is independent of the 
system size: e ~ 1//3. However g-breathers and normal forms predict e ~ 1/((3N 2 ) in that case. Thus derealization 
happens for much smaller energy densities than predicted by the IC criterion. 

The reason for the above discrepancy is the chosen resonance. Indeed, as we know from the perturbation theory of 
g-breathers, and from the analysis of normal forms, the leading order resonances are those between the third harmonic 
of the seed normal mode and the normal mode with three times smaller wavelength: A$h — w 3fc — 3wfe . Within 
the context of resonance overlap this was first recognized by Shepelyansky For long wavelength it follows from 
©, that we may approximate this frequency by A$h ~ &o- Comparing that with 6u> and using Uk « ko we find the 
correct threshold condition 



(49) 
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This condition is in excellent agreement with the scaling properties of strongly localized g-breathers, and with the 
scaling properties of resonant normal forms. We note, that the slightly more cumbersome analysis of the a-FPU case 
yields exactly the same conclusions (compare [f| and 0). 

VII. CONCLUSIONS 

We find that all three approaches - resonant overlap, periodic orbits (for strong localization) and resonant normal 
forms - yield identical scaling relations, if the proper resonance condition is identified and exploited. These scaling 
relations are depending on intensive quantities only. Consequently the above formulated FPU problems can be ex- 
plained in the following way. There are two time scales, since the excitation of a normal mode below the derealization 
threshold is close to an exact periodic orbit - a g-breather - which localizes exponentially in normal mode space. It 
will need a short time t\ to approach the g-breather. The excited trajectory will then spend a long time T2 in a 
neighbourhood of the periodic orbit. When approaching the derealization threshold, the two time scales merge, and 
equipartition is established on a single time scale. The thresholds depend on intensive quantities only, and thus the 
FPU effect is robust when considering the limit of infinitely large systems - but one has to remember to use intensive 
parameters. E.g. the often used scheme of exciting the mode with 50 = 1 but increasing the system size would lead to 
a change of the intensive quantity of a wave number fc or the frequency lu^ . An open question concerns whether the 
different scaling k c i as compared to k c i obtained for g-breathers, is also observed within the normal form approach. 
Another open point is the analytical estimate of the times t\ and and their relation to t c . 

Finally the generalization to higher spatial dimensions does not seem to bring about qualitatively new features, at 
least from the point of view of g-breathers. However, we think that much more detailed studies have to be conducted 
to be on the more safe side with such a statement (see [3l[ for some recent studies). 

To end with, we point out that Lichtenberg [33] performed recently an estimate of the (weak) stochasticity threshold 
for the a-model, similar to previous studies for the /3-model [24j . Surprisingly the results for the a-model are in full 
agreement with the scaling laws following from the analysis of g-breathers and normal forms. 

Acknowledgments 

A. P. thanks G. Benettin, L. Galgani and A. I. Neishtadt for useful comments and discussions on the subject. 



[1] E. Fermi, J. Pasta and S. Ulam, in Enrico Fermi - Note e Memorie (collected papers), vol. II (USA 1939-1954), Accademia 

Nazionale dei Lincei, Roma (Italy), and The University of Chicago Press, Chicago (USA) 1965, p. 978. 
[2] J. Ford, Phys. Rep. 213 271-310 (1992). 

[3] L. Galgani and A. Scotti, Phys. Rev. Lett. 28, 1173-1176 (1972). 
[4] J. Ford, J. Math. Phys. 2 387-393 (1961). 

[5] F. M. Izrailev and B. V. Chirikov, Sov. Phys. Dokl. 11, 30-32 (1966). 
[6] D. L. Shepelyansky, Nonlinearity 10, 1331-1338 (1997). 

[7] S. Flach, M. V. Ivanchenko, and O. I. Kanakov, Phys. Rev. Lett. 95, 064102/1-4 (2005). 

[8] S. Flach, M. V. Ivanchenko, and O. I. Kanakov, Phys. Rev. E 73, 036618/1-14 (2006). 

[9] M. V. Ivanchenko, O. I. Kanakov, K. G. Mishagin and S. Flach, Phys. Rev. Lett. 97, 025505/1-4 (2006). 
[10] O. I. Kanakov, S. Flach, M. V. Ivanchenko and K. G. Mishagin, Phys. Lett. A 365, 416-420 (2007). 
[11] T. Penati and S. Flach, Chaos 17, 023102/1-16 (2007). 

[12] S. Flach, O. I. Kanakov, K. G. Mishagin and M. V. Ivanchenko, Int. J. Mod. Phys. B, in press. 
[13] A. Ponno, Europhys. Lett. 64 606-612 (2003). 

[14] A. Ponno, in the Proceedings of the Cargese Summers School 2003 on Chaotic Dynamics and Transport in Classical and 
Quantum Systems, edited by P. Collet et al., Kluwer Achademic Publishers, Netherlands, 2005, p. 431. 

[15] A. Ponno and D. Bambusi, in the Proceedings of the International Conference on Symmetry and Perturbation Theory 
2004, edited by G. Gaeta et al., World Scientific Publishing, Singapore, 2005, p. 263. 

[16] A. Ponno and D. Bambusi, Chaos 15 015107 (2005). 

[17] D. Bambusi and A. Ponno, Comm. Math. Phys. 264, 539-561 (2006). 

[18] D. Bambusi and A. Ponno, in the Review Articles on Fermi-Pasta-Ulam Experiment grown out of the meeting FPU fifty 
years later held in Rome (Italy) on May, 7th 2004 (to appear on a Springer Book edited by G. Gallavotti); available on 
the web site |http ://ipparc o.romal.infn. it/index. html| 

[19] N. J. Zabusky and M. D. Kruskal, Phys. Rev. Lett. 15 240-243 (1965). 

[20] L. Berchialla, A. Giorgilli and S. Paleari, Phys. Lett. A 321, 167-172 (2004). 



14 



[21] F. Fucito et ai, J. Physique 43, 707-713 (1982). 

[22] J. Ford and J. Waters, J. Math. Phys. 4, 1293-1306 (1963). 

[23] A. Lyapunov, Probleme general de la stabilite du mouvement, published by the Mathematical Society of Kharkov, 1892 
(Russian); available in Editions Jacques Gabay, Paris (France) 1988 (French), reprint of the French translation from 
Russian appeared on the Annales de la Faculte des Sciences de Toulouse 2eme serie, tome IX, 1907. 

[24] J. De Luca, A. Lichtenberg and M. A. Lieberman, Chaos 5, 283-297 (1995). 

[25] L. Berchialla, L. Galgani, and A. Giorgilli, DCDS-A 11, 855-866 (2004). 

[26] J. A. Biello, P. R. Kramer, Y. V. Lvov, Stages of energy transfer in the FPU model, DCDS Supplements, 2003; spe- 
cial number devoted to the Proceedings of the Fourth International Conference on Dynamical Systems and Differential 
Equations, May 24-27, 2002, Wilmington, NC, USA, p. 113. 

[27] G. Benettin, R. Livi and A. Ponno, work in progress. 

[28] A. Giorgilli, S. Paleari and T. Penati, DCDS-B 5, 991-1004 (2005). 

[29] B. V. Chirikov, Internal Report n. 267 of the Nuclear Institute of the Siberian Section of the USSR Academy of Sciences, 

Novosibirsk, 1969; CERN Translation n. 70-41, Geneva 1971. 
[30] A. H. Nayfeh, Introduction to Perturbation Techniques (Wiley, New York, 1993). 
[31] G. Benettin, Chaos 15, 015108 (2005). 
[32] A. Lichtenberg, private communication. 



